SEMINARIO DI ANALISI MATEMATICA 
DIPARTIMENTO DI MATEMATICA 
DELL'UNIVERSITA' DI BOLOGNA 

Anno Accademico 1996-97 


Davide Guidetti 


SCHEMI ALLE DIFFERENZE FINITE E REGOLARITÀ 
MASSIMALE PER EQUAZIONI PARABOLICHE 
IN SPAZI DI FUNZIONI HÒLDERIANE 


22 maggio 1997 


Tecnoprint - Bologna 1997 


138 


Riassunto. Sono esposte alcune versioni per equazioni alle differenze 
finite di risultati di regolarità massimale per problemi parabolici in spazi 
di funzioni hélderiane. Vengono considerati gli schemi di Eulero e di 
Crank-Nicolson. Infine viene ottenuto un teorema di convergenza. 


Summary. Some finite difference analogs of results of optimal regularity 
for linear parabolic. boundary value problems in spaces of Hélder 
continuous functions are given; the backward Euler scheme and the 
Crank-Nicolson scheme are examined. 

We conclude with a result of convergence. 
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Scopo di questo seminario é presentare alcune versioni per equazioni alle dif- 
ferenze finite di risultati di regolarità massimale valevoli per problemi misti di tipo 
parabolico, ottenuti all’autore in collaborazione con S. Piskarev dell’Università di 
Stato di Mosca. A mia conoscenza la letteratura in proposito non é molto vasta. 
Potremmo dire che praticamente quasi tutto ciò che esiste già sull'argomento é 
opera della scuola di Sobolevskii a Voronezh. Segnalo a questo proposito la mono- 
grafia [1], che contiene essenzialmente applicazioni al solo caso di problemi su tutto 
R” e tratta la discretizzazione di problemi in cui la parte ellittica dell’operatore 
parabolico genera un semigruppo analitico nello spazio di base. É ben noto che cié 
non accade per il problema di Dirichlet se si lavora in spazi di funzioni h6lderiane 
ed é proprio questo caso che vorrei trattare. Mi limiterò a considerare il più sem- 
plice problema parabolico possibile, anche se molto di cié che diré potrebbe essere 
largamente generalizzato. Per questo faccio riferimento al preprint [2]. Premetterò 
ai risultati a cui accennavo una breve introduzione al metodo delle differenze finite 
nell’ambito delle equazioni paraboliche. 

Considererò dunque solo il seguente problema: 


Gt(t,a) = È4(t,2) + f(t,2), 
O0<t<T,0<x<L, 
u(t,0)=u(t,L)=0, 0<#<7T, 
u(0,x)=0, 0<xz<L. 


con T' e L numeri reali positivi. 
Siano M e N interi positivi; poniamo 7 := Th = L. Per n € {0,...,.M — 1},j € 
{0, ..., N}, una approssimazione di du(n7, jh) é, ad esempio, 


u((n+1)7,jh)- u(nr,jh) 
T ’ 


una approssimazione di 24 (n7,jh) per n E {0;...,M},j € {1,..,N- 1},é, ad 
esempio, 
u(nT, (jf+1)h)-2u(n7,jh)+u(n7, (jd —1)h) 
h2 i 
Poniamo ora, per 0<n<M-1,1<j<N-1, 
FP := f(n7,jh). 
In base alle considerazioni precedenti, consideriamo il sistema 
n+1 n n L_ n n . 
Sil VU pr, per0<n<M-1,1 <j<n-]1, 
UT =U%r=0per1<n<M, (2) 
U? =0 per0<j< N. 


Introduciamo ora le seguenti notazioni: sia 7 un insieme ed E uno spazio di Ba- 
nach con norma ||.||; indichiamo con B(7, E) l'insieme delle applicazioni da T a E 
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limitate. Omettiamo di indicare E se E = C. Se F € B(T, E) e j € T scriviamo 
F; o Fi invece di F(j). Infine, se G € B(Z, E) poniamo ‘ 


||G]|3(,5) = sup ||G;||. 
jEr 


Allora, evidentemente, U” é un elemento di B (Cy); con Tn = {L....N- 1}. Per 
V e B(In), j € {1,...N — 1}, poniamo 


Vizi 2VG+ Vi 


AnVj = TO 


(3) 
con la convenzione che V; = 0 se j € {0,N}. Possiamo allora riscrivere (2) nella 
forma casi 
uf =U" = A,U"+F",n=0,..,M-1, (4) 
U©=0 


che può essere evidentemente risolto ricorsivamente. 

Ora, una proprietà fondamentale che deve essere posseduta da uno schema di dis- 
cretizzazione 6 la stabilità: nel caso in questione, questa proprietà può essere for- 
mulata nel seguente modo: esiste C > 0 indipendente da 7 e h tale che per ogni 
FP=a(P3,. EV 4) 


max ||U"||sam < C_max_ [|F"]l8cwm- (5) 


0<n<M O0<n<M-1 


É ben noto che lo schema (3) non é incondizionatamente stabile, nel senso che la 
stabilità viene a cadere se non si presuppone qualche relazione tra 7 e h. Una 
relazione che garantisce la stabilità é la maggiorazione 


1 


Per questo e altri risultati di tale genere, si veda il lungo articolo [3]. 
Naturalmente, (4) (detto schema di Eulero esplicito) non é l’unico schema di dis- 
èretizzazione possibile per (1). Consideriamo il cosi detto "schema di Eulero im- 
plicito” o schema di Rothe: 


U-UP2 — AU" + F",n=1,...M, (7) 
U®=0 

Osserviamo che il calcolo della soluzione di (7) é molto pit laborioso del calcolo 

della soluzione di (4), in quanto richiede l'inversione di una trasformazione lineare: 

si vede abbastanza facilmente che per ogni 7 e ” l'operatore lineare 1 — TÀ, é 

invertibile in B(Zy) e per la soluzione si ha, dato n € {1,...,.M}, 


UU" SN hl * FP (8) 


i=1 
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La maggiore difficoltà di calcolo della soluzione di (7) é tuttavia compensata dalla 
proprietà di stabilità incondizionata dello schema. 

Un altro schema interessante é il così detto schema di Crank-N icolson, che é a metà 
strada tra i due schemi di Eulero: 


UPZUTT — 1(ALU® + A,U"-1)4 pon 1,...,.M, (9) 
U.=0. 


In questo caso la soluzione di (9) é data da 


n 


U" = 7(1- ZAN) YO r(rA,) A, (10) 


l=1 


ove r(2)=(1+5)(1- 3). 
L’interesse dello schema (9) sta nel fatto che é pi accurato dei due precedenti nel 
seguente senso: sia u la soluzione di classe C® su [0,7] x [0, L] di (1) (ovviamente 


supponiamo che (1) ammetta una soluzione con questa regolarità). Poniamo, per 
n € {0,...,M},j € 40; IV 


uj := u(nr,jh) (11.) 
Allora, si ha, per lo schema (7), 


u — ya! 
—_ = Anu"+F" + O(7+ h2), (12) 


T 


ponendo F? = f(n7,jh). Invece, per lo schema (9), 


u"-u"1 1 

2 = gAn(u"+u"7)+ F"+0O(72+h?), (13) 
ove Fi = f((n- 1)r, ih). Applicando la formula di Taylor, nel caso in cui si abbia 
stabilità, la maggiore accuratezza di una schema Barantisce, per piccoli valori di 
T e h, una migliore approssimazione della soluzione di (1) mediante la soluzione 
dello schema alle differenze. É ben nota la stabilità incondizionata dello schema 
(9). Rilevo comunque il fatto per un generico problema di evoluzione, la stabilità 
dello schema di Crank-Nicolson 6 un problema difficile. Ricordo solo il seguente 
recente risultato, conseguenza quasi immediata di stime dovute a vari autori (vedi 
[4] e (5): 


Teorema 1. Sia E uno spazio di Banach, A un operatore lineare in E non neces- 
sariamente limitato, con il risolvente contenente {z € C — {0}||Arg(z)] < w}, per 
qualche w €]Z,t[. Valga inoltre, per 2 nel settore detto, 


II(z- A) lc < Ma, (14) 


con M > 0. Consideriamo il problema 


T 


Wper> _ i n n-1l mM pu 
—7— = 3(AU" + AU") + F°,n=1,..,M, (15) 
UO -0, 
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con F*,...,F" € E. Allora, (15) é univocamente risolubile ed esiste C > 0 dipen- 
dente solo da w.M.T tale che . 


n < é D_ 
pax "Ile SC mex IP"ile (16) 


E’ chiaro che se si discretizza il problema anche nello spazio (come nel caso 
particolare che ci stiamo limitando a considerare), la validità uniforme nel passo di 
discretizzazione spaziale delle ipotesi del teorema 1 implica la stabilità. 

Passiamo ora a stimare norme pit forti della soluzione U dello schema alle differenze; 
introduciamo a tale scopo un po’ di notazioni. 

Sia T C Z e sia V € B(Z, E) dove E é uno spazio di Banach. Sia poi @ €]0, 1[. 
Poniamo 


Mora) = sup [(k — |a — Vill. (17) 
j,kET,j<k 
IVilcea,a = max{I|Vilec.a): Moga.a)} (18) 


Introduciamo ora la differenza divisa in avanti: poniamo Ti := {j € T|j+1 € 2} e, 
per j € T,, h> 0, V € B(Z,E). 


o, = et. 


(19) 
Pit in generale, sia per / intero positivo, Z, := {j € T]j+1,....j+! € T}. Perj € T 
é ben definita la differenza divisa di ordine / 


l 
S: l a 
d,V;:=h') (en pri (20) 
i=0 
Per V € B(I, E), m intero positivo, poniamo 


. _ i 
|Vllce (e) = RA [o V]lBa,5) (21) 


e, se m é un intero non negativo e 8 €]0, 1[, 
IVilorto o = maxf{IVilop aa): (0h Vor r.a)} (22) 


Vediamo ora un risultato di regolarità massimale parabolico valevole per il problema 
(1); preliminarmente, definiamo soluzione stretta di (1) una funzione di dominio 
[0, 7] x [0, L] continua assieme a qu du, DI, Si ha allora: 

Teorema 2. Consideriamo il problema (1) e sia 8 €]0,2(-{1}; le seguenti con- 
dizioni sono necessarie e sufficienti affinché (1) possieda una soluzione stretta tale 
che t + u(t,.) € B([0,T]; C2+°([0, L})),t > $t(t,.) € B((0.T]; C°({0.L]))): 

(a) f e C(,T} xD); 

(6) #> f(t..) € B(0,T}: 0°((0, L)); 

(c)t > f(t,0) et— f(t.L) appartengono a C([0. T)); 
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(d) f(0,0) = f(0,L) =0. 


Vogliamo dare qualche versione del teorema 2 valevole per schemi alle differenze. 
Cominciamo dallo schema di Eulero implicito. Introduciamo alcune notazioni: per 
M intero non negativo poniamo Ty := {0,1,...N-1,N} = ZyU {0, N}. Se 
V € B(Iy) indichiamo con V l’elemento di B(Ty cosi definito: 


v=f Viif1<j<N-1, 


0 if je {0, N}. ini 


Useremo spesso la notazione £V in alternativa a V. Se G € B(In), indichiamo con 
RG la sua restrizione a B(Zy). 

Il seguente risultato preliminare può chiarire il significato delle ipotesi del risultato 
principale che seguirà: 


Lemma 1. Sia Ax l'operatore introdotto în (3) su B(In): per V € B(Iy), 
AnVj = h72(Vi-1- 20; + Via). 


Allora: 

a) esiste 6 > 0 indipendente da h tale che {: € C- {0}j]Arg(A)] < t}U {A e 
b) per ogni $ € [0,r[ esiste M(é) > 0 dipendente solo da $ tale che, se A #0 e 
|Arg(A)|< @ oppure |A| < 6 


IA — 44) [leeam) £ M(O(1 +14); (24) 


c) se 0 €]0,2(-{1}, [él < 7, A e C- {0}, |Arg(A)] < @ 0 |A < 6, per ogni 
GE B(INn) 


IATTIE(A — 414)! AGI|og em + IE — 44) RG] |c2+0 7, 


< (0, 9)[IG]lcey) + |A? max{[Gol, |Gw]}}, (25) 


ove c(9, $) dipende solo da 0 e $; 
d) esistono due costanti positive c1(0) e c2(9) dipendenti solo da 0 tali che per ogni 
Fe B(In) 


(0Feem.pam)g £ IEFilogam £ 2(0F een DA), (26) 
Per F € B(Iyw) abbiamo posto 


Foa») = max{||Fllaw), |\An-Fll8(w) (27) 


e ||Fl am), D(A4,) a __ indica la norma nello spazio di interpolazione reale. Natu- 
00 


ralmente, per À fissato, queste due norme sono equivalenti alla norma- Il.llacy) € 
(26) stabilisce l'equivalenza uniforme tra due classi di norme. 

Veniamo ora al risultato principale relativo allo schema di Eulero implicita: nella 
versione che dé é un caso particolare di 4.3 in [2]: 
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Teorema 3. Sia B, un operatore lineare in B(Iyn) soddisfacente le seguenti con- 
dizioni, per certi R > 0,9 €]T.m[, M > 0,0 €]0,2[-{1}: ‘ 
(a) {A € CI|A| > R. |Arg(A)| < do} C p(Ba) e, per A in questo insieme, 


IA A) lean) £ MA + (28;) 
(b) per ogni F € B(In) 


EF|lceam £ MIF|lx,D(B))g (29) 


(c) per ogni G € B(Ty),A € C con |A > Re | Arg(A)| £ do 


(+ DNIEO — ARG oper) + IE — AVRO cen) 


< MI|{Glloeew) + (1+1A))? max{[Gol, [Gw]}}. (30) 


Si consideri il problema 

un-u?"! _ B,U" + RG",n=1,..,M, 

5 (31) 
u'=0. 


con G € B({0,1,..., M}, B(Zy) tale che G°=0 
Allora, se TR < 1, (31) é univocamente risolubile; se U € B({0,1,..., N}, B(Zn)) 
risolve (31), per n = 0,1,...,.M 


|EU"]|cz+o gw £ Aa, IG" |lce Em) 


(ma — na )7) 8168? — GEN + _ ax (ma —ma)7) FIG - GR 
(32) 


ove € é una costante positiva dipendente soltanto da 9, R, do,M,T e indipendente 
dar, heG. i 


Osserviamo che in base al lemma 1 gli operatori A, verificano le condizioni del 
teorema 3 con costanti indipendenti da h. Vale perciò la conclusione del teorema 3 se 
si scelgono G!,...., GM in maniera tale che il secondo membro della disuguaglianza in 
(32) resti limitato da una costante indipendente da 7 e h. Se il dato f nel problema 
(1) soddisfa le ipotesi del teorema 2, é facile vedere che, se poniamo, fissati 7 e h 


È max 
0<n1<n2<N 


Gj = f(n7,jh) (33) 
per 1<n<Me0<j<N,ilsecondo membro di (30) é minore o uguale di 


: se) -|f(.0)] ..L è 
A INF: ilcoro.ap + (il gr + IE D Ig po BI 
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Un altro operatore alle differenze su B(Zw) che soddisfa le condizioni del teorema 
3 uniformemente rispetto al parametro h, purché si scelga il parametro p sufficien- 
temente grande, é il seguente: per U € B(In), 1<j<N-1, poniamo 


-— 00,2 + 1920;_, = 1020, se 2 < j < IUo= 2, 
Choi = 0?Uo — 5 [02U, — 202U,41+ 020,42] sej= 1, 
RÙ,.-5 _ D(02Un-p-4 _ 2000, -pa + RÙ,,-- pal sej=n_- 1. 
(35) 
Tale operatore rappresenta una discretizzazione pit accurata della derivata seconda 
rispetto a x rispetto ad A, valendo per u € C*([0, L]), 


u"(jh) = Apu? + 0(h2), (36) 


u"(jh) = Ch, pu? + 0(h8). ° (37) 


Passiamo ora allo schema di Crank-Nicholson. Il risultato corrispondente al teorema 
3 é il seguente teorema 4 (caso particolare di 5.4 in [2]): 


Teorema 4. Supponiamo che valgano tutte le ipotesi del teorema he e, in più, 
d) ||rBxllcxaaw) £ S, con S>0; 
e) se |A] > 25 


EA — TA) RG] cory £ MIIG]lcg ay +77’ max{Go, Gw}] (38) 


per ogni G € B(In). 
Si consideri il problema 


T 


U2U” _ 1(ByU" + BU"-1) + RG",n=1,..,M, (39) 
U®=0.. 


Allora, se 2rR < 1, (39) é univocamente risolubile e vale la conclusione del teorema 
3. 


Qualche parola di commento su questo risultato: la condizione discriminante é 
d); essa impone nella pratica una relazione tra il passo spaziale e il passo temporale. 
Infatti, consideriamo, ad esempio, l’operatore Ax; si ha 


Akllecsawy = 487. (40) 


d) impone allora 
T<ah?, (41) 


si può verificare che il teorema 4 fornisce stime di stabilità sia nel caso degli operatori 
Ah, sia nel caso degli operatori B, purché si scelga p sufficientemente grande. In 
questo caso il risultato é stato provato solo nel caso @ €]0, 1[. 

Veniamo ora al problema della convergenza delle soluzioni dei problemi alle dif- 
ferenze alla soluzione del problema differenziale; cominciamo con lo schema di Eu- 
lero implicito; solitamente si presuppone una certa regolarità della soluzione. 
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Teorema 5. Supponiamo che (1) ammetta una soluzione stretta u cont + ult, .) € 
C1+8((0,T]; C([0, L])) n B([0.T}; C+9({0, L})), t + 24(t,.) e B((0,7]: C°([0,1))), 
con 8 €]0, 1], in]0, 2[; sia U la soluzione del sistema (29), con Br = An e G7 = 
f(nr,jh); allora, per ogni £ € [0, 0], e > 0 esiste c > 0 dipendente solo da u, £, e tale 
che, per n = 1,...,.M 


lu — EU" ]ogrea S AP + 09)IT$7<, (42) 


u — un! Ut — un! 


T T 


Stime dello stesso genere (vedi [2], 6) valgono per lo schema di Crank-Nicolson. 
ponendo in tal caso By = Ck p; con p > po, prendendo 7 < ah? per qualche a > 0. 
In questo caso, tenendo conto di (11), é conveniente porre FE = f((n- 3)7, 3h). Se 
la soluzione di (1) é sufficientemente regolare, per certi valori di £ (vedi il teorema 
5) le soluzioni dei problemi discretizzati convergono più rapidamente alla soluzione 
di (1). l 
Concludiamo osservando che risultati dello stesso genere si possono ottenere per 
problemi misti di Cauchy-Dirichlet arbitrari in dimensione uno e anche per certi 
problemi. in dimensione superiore (vedi [2]). 


llesam £ ce(r8f + n9y-i—e. (43) 
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